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THE  ELASTIC  MODULI  OF  HETEROGENEOUS  MATERI/  LS 

by 

Zvi  Hashin"*' 

1.  INTRODUCTION 

The  present  work  is  concerned  with  the  determination,  by  theoretical 
analysis,  of  the  elastic  moduli  of  a  randomly  heterogeneous  material. 

It  is  assumed  that  the  material  may  he  adequately  described  by  an 
elastic,  homogeneous  and  isotropic  matrix  the  moduli  of  which  are  known, 
in  which  particles  of  another  elastic  homogeneous  and  Isotropic  material 
are  imbedded,  the  moduli  of  which  are  also  given.  Assuming  furthermore 
that  the  volume  concentration  of  particles  is  uniform  and  that  the  material 
may  accordingly  be  regarded  as  quasi-homogeneous  the  problem  is  to  find 
expressions  for  the  effective  elastic  moduli  of  this  heterogeneous  material. 

Numerous  papers  on  the  determination  of  bulk  properties  of 
heterogeneous  materials  have  been  published.  The  first  of  these  was  by 
Einstein  (1906,  1911)  in  which  the  viscosity  of  a  suspension  was  determined, 
assuming  that  it  may  be  described  by  rigid  spheres  suspended  in  a  viscous 
fluid  and  that  the  volume  concentration  is  so  small  that  particles  do  not 
interact.  The  case  of  dilute  concentration,  assuming  that  the  particles  are 
spherical,  has  been  solved  for  a  varied  of  materials:  Liquid  droplets  in 
another  liquid,  Taylor  (1932);  Elastic  particles  in  viscous  fluid;  Froehlich 
and  Sack  (1946);  Empty  holes  in  elastic  solid,  Mackenzie  (1950);  Rigid 
particles  in  elastic  solid,  Hashin  (1955);  Elastic  particles  in  another  elastic 
material,  Eshelby  (1957),  Hashin  (1958);  Viscous  liquid  inclusions  in  elastic 
solid,  Oldroyd  (1956).  Metiiods  of  analysis  for  small  concentration  have 
been  recently  applied  b>  Budiansky,  Hashin  and  Sanders  {I960)  to  the 
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theoretical  determination  of  the  early  plastic  behavior  of  polycrystalline 
materials. 

In  the  case  of  small  volume  concentration  the  fimdamental  assumption 
is  that  partioles  do  not  interact.  The  effective  moduli  can  then  be  determined 
in  a  way  outlined  below  (Section  2), 

The  problem  of  finite  concentration  is  much  more  complicated.  There 
is  at  present  little  hope  of  rigorously  solving  the  problem  of  a  medium  in 
which  there  are  many  interacting  inclusions.  The  difficulty  of  this  kind  of 
problem  is  well  illustrated  fay  Sternberg’ s  and  Sadowslq^'  s  (1952)  solution  of 
the  axisymmetric  problem  of  the  theory  of  elasticity  containing  only  two 
spherical  cavities.  [See  also  Miyamoto  (1956).] 

Most  of  the  work  which  has  been  done  up  to  now  in  the  field  of  finite 
concentration  has  been  concerned  with  the  extension  of  Einstein' s  formula, 
previously  mentioned,  to  higher  concentrations.  This  has  resulted  in  a  large 
number  of  different  formulae  which  are  sometimes  in  direct  contradiction. 

Up  to  now  no  formula  has  been  derived  which  is  theoretically  well  fotmded 
and  fits  experimental  results  for  the  whole  range  of  volume  concentration. 

A  review  of  these  investigations  is  not  within  the  scope  of  this  paper. 
Summaries  of  these  may  be  found  in  Frisch  and  Simha  (1956)  and  Reiner  (1958). 

Instead  of  an  attempt  to  find  an  exact  formula  for  the  moduli  of  a 
heterogeneous  body  this  work  will  be  concerned  with  the  construction  of 
approximate  upper  and  lower  bounds  for  the  moduli.  Obviously  such  bounds  are 
of  practical  value  only  if  tney  are  close  together.  It  will  be  found  that  the 
expressions  for  thfi  bulk  modulus  bounds  coincide.  In  the  case  of  the  shear 
modulus  they  are  mostly  close  together,  as  is  illustrated  by  a  specific  case 
which  has  been  solved  numerically  in  this  paper.  Also  a  simple  expression 


is  derived,  which  is  always  smaller  than  the  upper  bound  and  larger  than 
the  lower  one  and  can  thus  be  used  as  a  good  approximation  to  the  shear 
modulus  whep<*yer  the  bounds  are  close. 

Apparently  tiie  first  bounds  on  an  elastic  modulus  of  a  heterogeneous 
material  have  been  given  In  a  recent  paper  by  Paul  (1959).  Bounds  for  the 
Young’ s  modulus  of  a  heterogeneous  material  were  obtained  by  using  the 
variational  theorems  of  the  theory  of  elasticity  and  taking  as  an  admissible 
stress  system  (see  Section  2)  the  same  simple  tension  in  matrix  and  particles 
and  a  simple  tension  deformation  for  an  admissible  displacement  field.  While 
these  bounds  have  the  advantage  of  being  exact  they  are  far  apart  because  of 
the  simple  admissible  stress  and  displacement  fields  chosen.  Thus  for 
example  for  the  same  material,  for  which  a  numerical  solution  has  been 
given  in  this  paper,  [  experimental  results  for  this  case  have  been  obtained 
by  Nishimatsu  and  Gurland  (1959)  and  will  be  given  below]  Paul  found  that 
for  a  concentration  of  50  per  cent,  the  upper  bound  for  Young’ s  modulus 
was  43  per  cent  higher  than  the  lower  bound. 

In  the  following  general  expressions  for  the  elastic  moduli  of  a  quasi- 
homogeneous  heterogeneous  material  will  be  developed,  involving  only  the 
stresses  or  strains  inside  the  particles.  This  is  done  by  considering  the 
change  in  strain  energy  in  a  loaded  homogeneous  body  due  to  the  insertion 
of  non-homogeneities. 


bounds  approximately  two  geometrical  approximations  are  made.  It  is 
assumed,  as  in  small  concentration  theory,  that  the  particles  are  spherical 
and  moreover  that  the  action  of  the  whole  heterogeneous  material  on  ai^r  one 
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particle  is  transmitted  through  a  spherical  shell  which  lies  wholly  in  the 
matrix. 

It  is  believed  that  this  approximation  is  a  close  one  and  comparison 
of  theoretical  and  experimental  results  strengthens  this  belief. 


2.  GENERAL  THEORY 

Let  an  elastic  homogeneous  and  isotropic  body  of  arbitrary  shape 
(Fig.  la)  be  subjected  to  <:vc*face  tractions  which  are  associated  with  a 
homogeneous  stress  tensor 


where 


Jo)  Jo) 

I’i  =  o-y  nj 


n<0) 


(2.1) 


i;  j  =  1,  2,  3  (a  repeated  subscript 
denotes  summation) 


*i 

(o) 


are  the  components  of  the  surface  stress  vector 

constant  stresses 
Uj  components  of  the  normal  to  the  outer  surface 

The  strain  energy  Uq  stored  in  the  bo^  is  then  given  by  either  of 


the  formulae 
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where  the  cT..  and  the  strains  e,.  are  connected  by  Hooke’s  law 
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If  part  of  the  body  is  replaced  by  an  inclusion  of  another  material, 
the  elas^tic  moduli  of  which  are  Kp,  Xp,  Gp  (Fig.  lb)  and  the  same  surface 
tractions  are  applied,  then  it  has  been  shown  by  Eshelby  (1951,  1956)  tha.t 
the  difference  in  elastic  energy  stored  in  the  two  cases  is  given  b7 

^  uH  =  =  I  \_  (Tj®^  -  Tj  u[°^)  dS  (2.10) 

^(S) 

where  S  is  the  surface  of  the  includon,  fj  and  u^  are  the  stress  and 
displacement  vectors  on  the  surface  c£  the  inclusion. 

Formula  (2.10)  is  general;  it  is  valid  for  any  inclusion  shape  and  any 
boundary  tractions. 

An  alternative  useful  form  for  (2.10)  has  been  given  by  Eshelby,  this  is; 


^  =  2  f *  V  ^  +  2  (G^  -Gp)  e[°^  ey]  dV  (2.11) 

This  may  be  rewritten  in  terms  of  stresses  by  using  (2.7) 

1  f  .  °m-Op  (0)  1 

2  9K,nKp  “« 


where  V  is  the  volume  of  the  inclusion.  Again  (2.11)  is  a  special  case  of 
Eshelby' s  result  which  is  valid  for  a  nonisotropic  and  nonhomogene ous 
inclusion.  When  instead  of  surface  tractions,  surface  displacements  are 
prescribed  formulae  (2,10),  (2.11)  and  (2.12)  change  in  sign.  If  when 
displacements  are  prescribed,  is  denoted  by  S  then  Eshelby  has  shown 
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2  J-  '’^i  "i  ^i^i  )  dS 

(S) 


(2.13) 
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The  derivation  of  (2.10),  (2.11)  and  (2.12)  remains  unchanged  when 
instead  of  one  inclusion  an  arbitrary  number  N  of  them  is  introduced  bito 
the  elastic  body  and  surface  tractions  are  held  fixed.  (Fig.  Ic) 

Then  the  contribution  of  the  n^  inclusion  is: 


(2.14) 


(2.15) 


Analogous  results  hold  for  (2.11),  (2.12)  and  (2.13).  In  the  following, treatment 

of  the  case  where  tractions  are  prescribed  will  be  called  the  stress  approach. 

when  displacements  are  prescribed  •  the  displacement  approach. 

The  bulk  and  shear  modulus  of  an  elastic  body  containing  a  large 

nxunber  of  inclusions  of  another  material  will  now  be  given  in  terms  of  energy 

(o) 

expressions.  It  is  first  assumed  that  the  T^  are  equivalent  to  isotropic 
tension 
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Then  from  (2.12),  for  the  n“  inclusion 
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When  instead  of  isotropic  stress,  isotropic  radial  displacement 
(O)  6(0) 

“  ■'  g '  '  Xj  is  prescribed  on  the  boundary,  then  from  (2.10),  (2.11)  rx.d 


(2.13),  and  proceeding  in  an  analogous  way 
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dV 
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For  the  homogeneous  body  without  inclusions 
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(2.20) 


If  it  is  assumed  that  the  heterogeneous  material  is  quasi -homogeneous  with 
bulk  modulus 


2  9K* 


(2.21) 


where  (2.21)  may  be  regarded  as  a  definition  of  the  "effective"  modulus. 
Then, 


n=N 


(O’) 


n=x 

which  on  using  (2.20)  and  (2.21)  may  be  written  as: 


(2.22) 


1  _  1 


(2.23) 
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where  i  U'  '  is  given  by  (2.22)  and  (2.18). 

It  should  be  borns  in  mind  that  (2,23)  is  valid  only  when  the  actual 
are  used  in  (2.18). 

Using  now  the  displacement  approach  the  following  expressions 
are  obtained: 


1 
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U^’ '  1^'*°’''^  <2-2^> 
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U  =  —  K  6  '  y  \^*25; 

jt 

* 
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2 

Analogous  formulae  can  be  obtained  for  the  shear  modulus  G  .  It  is 

convenient  to  apply  the  stress  vector  field 

r. 

T  “  f'n-. 

Aj  Cli2 

(0) 

T2  =  tn^  (2.27) 

(0) 

Tj  =  0 

which  is  equivalent  to  the  state  of  stress, 

»=0  s‘?  =  s<°>=T 

<»)  _  (0)  ,(o)  (0)  (0) 

^li  S22  -  ^33  =  Sjj  -  ^Xi  -  0 

which  is  a  pure  shear. 

Then  from  (2.12); 
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Also, 

2  Gm 

=  -i-j  r^V  (2.31) 
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and  using  these  expressions  in  (2.22) 


G" 
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1 


(2.32) 


Proceeding  in  an  analogous  way  when  homogeneous  shear  displace¬ 
ments  are  prescribed  on  the  boundary,  given  by 
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the  following  results  are  obtained: 
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where 


u<‘>  =  1  G*  f  V 
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The  expressions  for  the  effective  moduli  depend  on  the  actual  stresses 
or  strains  in  the  inclusicms.  In  the  case  of  small  volume  concentration,  when 
it  is  assumed  that  there  is  no  interaction,  these  may  be  determined  by 
assuming  that  each  inclxision  is  in  an  infinite  medium  where  the  stresses  or 
strains  at  infinity  are  the  ones  applied  to  the  surface  of  the  body.  Solving 
the  boundary  value  problem,  the  expressions  for  the  case  of  small  concentration 


are  easily  obtained.  This  procedure  is  clearly  impossible  for  finite 
concentration.  On  the  other  hand  bounds  on  the  moduli  may  be  obtained 
using  the  variational  theorems  of  the  theory  of  elasticity. 

The  theorems  of  minimum  complementary  energy  and  of  mtnimiun 
potential  energy  will  be  here  used  in  their  following  special  forms: 

(a)  When  tractions  are  prescribed  over  the  entire  surface  of  an 
elastic  body  and  the  body  forces  vanish,  then  of  all  sets  of  stresses,  satisfying 
the  equilibrium  conditions  and  boundary  conditions,  the  actual  state  of  stress 
minimizes  the  strain  energy 

(b)  When  displacements  are  prescribed  over  the  entire  surface  of  an 

elastic  body  and  the  body  forces  vanish,  then  of  all  sets  of  displacements, 
satis^ing  the  boundary  conditions,  the  actual  displacements  minimize  the 
strain  energy  [  c.f.,  e.g.  Sokolnikoff  —  Mathematical  Theory  of 

Elasticity  (1956),]  Thus  for  the  special  cases  mentioned,  both  theorems 
reduce  to  a  minimum  principle  for  the  strain  energy- 

In  what  follows  a  field  of  stress  or  displacements  belonging  to  the 
sets  described  in  (a)  or  (b)  respectively  will  be  called  admissible. 

If  in  the  stress  approach  any  admissible  state  of  stress  is  chosen, 
the  strain  energy  associated  with  it  may  be  computed  from  (2.12)  and  (2.22). 

If  the  expression  S  thus  determined  is  denoted  by  6  and 
is  the  actual  strain  energy,  then  from  (2.22)  and  principle  (a): 

u(o')  ^  (2.38) 

0 

If  tow  in  the  displacement  approach  an  admissible  field  of  displace* 
ments  is  chosen,  the  strain  energy  is  computed  from  (2.11)  and  (2.22)  and 
the  expression  S  is  denoted  l?y  ^  -*  then  from  principle  (b) 
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u(^)  s  (2.39) 

The  bounds  on  the  bulk  modulus  are  obtained  in  the  following  way: 
An  admissible  state  of  stress  is  chosen,  so  that  on  the  boundary  S  of  the 
body  shown  in  Fig.  Ic  (2,16)  is  satisfied,  and  J  is  determined. 
From  (2.38),  (2.20)  and  (2.21) 


which  may  be  rewritten  in  the  form: 


(2.42) 
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HETEROGENEOUS  MATERIAL  DIVIDED  INTO 
COMPOSITE  ELEMENTS 
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where  Vj^  is  the  volume  of  the  inclusion  and  V  the  volume  of  the  whole 
body.  By  uniform  concentration  the  following  is  understood:  when  choosing 
an  arbitrary,  not  too  small,  volume  element  of  the  body  the  fractional  volume 
of  inclusions  in  it  is  expected  to  deviate  from  c  only  by  a  small  amount. 

A  system  of  tractions  derived  from  homogeneous  stresses  is 
now  applied  to  the  surface  of  the  body.  Each  inclusion  is  then  imagined  to 
be  surroimded  by  a  surface  Sn,  enclosing  a  volume  Vn,  so  that  for  every 
inclusion 


It  is  possible  to  construct  these  surfaces  Sq  in  an  infinity  of  ways  and  it 
will  be  assumed  that  this  has  been  done  so  as  to  approach  spherical  surfaces 
as  nearly  as  possible.  In  the  following  any  element,  containing  an  inclusion 
and  enclosed  by  will  be  called  -  a  composite  element.  (Fig.  3a) 

In  order  to  choose  a  state  of  stress  which  satisfies  equilibrium 
and  boundary  conditions  it  will  be  assumed  that  every  surface  Sjj  of  a 
composite  element  is  loaded  by  the  tractions  acting  on  S  .  If  then  the 
boundary  value  problems  for  the  composite  elements  are  solved,  the  stresses 
thus  determined  fulfil  the  necessary  conditions  for  an  admissible  stress 
system  .  In  order  to  solve  the  problem  it  is  assumed  that  the  composite 
element  may  be  approximated  by  two  concentric  spheres  so  that  volumes 
are  preserved  (Fig,  3b).  Such  an  element  will  be  called  -  a  spherical 
composite  element. 

For  the  displacement  ^proach  it  is  assumed  that  the  surface  displace¬ 
ments  of  S  are  linear  and  that  the  surface  displacements  of  each  Sq  are  of 
the  same  linear  form.  When  the  corresponding  boundary  value  problem  for 
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the  composite  element  is  solved  the  resulting  field  of  displacements  will 
be  admissible. 

3.  APPROXIMATE  BOUNDS  FOR  THE  BULK  MODULUS 


(a)  Lower  bound  by  stress  approach. 

The  spherical  composite  element  is  given  by  a  sphere  of  radius  b^ 


concentric  with  a  spherical  particle  of  radius  a^ .  By  definition 


^n  ^ 


(3.1) 


where 


(3.2) 


The  outer  surface  S  of  the  heterogeneous  body  is  loaded  by  a 


constant  radial  stress 


(3.3) 


Then  the  same  .stress  is  applied  to  the  spherical  surface  r  =  b .  (hi  the 
following  the  subscript  n  will  be  dropped.) 


The  elastic  moduli  of  the  inclusion  are  Gp .  Those  of  the 
spherical  shell  -  Km;  Gm* 

The  problem  is  one  of  radial  symmetry  the  general  solution  of 
which  is  [  (Love),  p.  142]: 


The  four  unknown  constants  may  easily  be  determined  from  the  four 
boundary  conditions: 
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From  the  preceding  analysis  it  follows  that  only  the  stresses  inside 
the  inclusion  are  of  interest.  It  is  found  that: 


(3.9) 


I/* 


A,- 


\  -i 


Then 


r  ~  ^  %  "S 


From  (2.18),  (3.1),  (3.2),  (3.9)  and  (3.10) 


=  I  . 


From  (2.40)  and  (2.47), 


4  Gjq  -f  3  Km 


(3.10) 


-"^m  (4Gm+3  Kp)  - 40^  (Km  "Kp)  c 

(3.11) 


Ki  = 


^  (Km  “  Kp)  (4  Gm+  3  Km)c 

""  Km  (4Gm  +  3Kp)  -4Gm(Km  “Kp)  c 


(3.12) 


(b)  Upper  bound  by  displacement  approach. 

It  is  here  assumed  that  the  boundary  displacement  of  S  is  purely 


radial 


.(o)  _  £ 


34  =  exi 


(3.13) 


Accordingly  the  displacement  applied  to  the  boundary  of  the  spherical 


composite  element  is 


(o) 

Up  =  er 


(3.14) 


The  boundary  value  problem  is  solved  exactly  as  before,  the  only 
difference  being  in  that  equation  (3.8)  has  to  be  replaced  l?y  (3.14)  for  r  =  b . 
It  is  then  found: 


'r 


V.  A' 


V. 

\  1. 


..  \ 
K" 


,4 


Ap=  t 


_ 4  Gm  3  Km _ _ 

4Gni  +  3Kp  +  3{Kni  -  Kp) 

D 


(3.15) 


Bp=  0 


(3.16) 


From  (2.19),  (3.1)  and  (3.15), 


“I' 


_ 4  Gm'*'  3  Xm 

4Gj„  +  3Kp  +  3(Kiji  -  Kp)  c 


and  from  (2.42)  and  (2.47) 


K3  =Kin  +  (Kp  -  Km) 


(4Gm  +  3Km)  c 
4Gm  +  3Kp  +  3(Km  -Kp)c 


(3.17) 


(3.18) 


Comparison  of  K^  given  by  (3.12)  and  Kf  given  by  (3.18)  shows  that  the 
upper  and  lower  bounds  coincide. 


k:  =  Kt  =  K^ 


(3.19) 


As  the  bounds  are  appro.ximate,  this  does  not  necessarily  show  that  the 
expression  for  R  is  exact. 

The  expression  (3.18)  can  be  rewritten  in  the  form 


=  1  + 


m  '  “m 


It  can  also  be  rearranged  into  the  S3nnmetric  form 


1  -  ^ 


(3.20) 


(3.21) 
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where  f  is  given  by 

,  4Gr 


.  2(1  -  2  Vm) 


1  + V* 


m 


and  Vjjj  is  the  Poisson's  ratio  of  the  matrix.  Generalization  to  the  case 
when  the  particles  are  of  different  kinds  is  obvious.  Let  there  be  k  kinds 
of  pai’ticles  imbedded  in  the  matrix.  Let  be  the  bulk  modulus  of  the  i*^ 
particle  kind  and  c^  its  voliune  concentration.  Then  (3.20)  is  generalized 
to: 


i=k 


1+3(1-V^) 


.,(i) 


K 


-  1  c 


m 


2(1-2V„,)  +  (1+V^m) 


L^^m 


(3.22) 

where 

i=k 

c  =  ^  Ci 
i=l 

4.  APPROXIMATE  BOLTroS  FOR  THE  SHEAR  MODULUS 

The  method  is  essentiaUy  the  same  as  for  the  bulk  modiilus.  A 
homogeneous  shear  stress  or  homogeneous  shear  displacement  is  applied 
to  the  boundary  S  and  also  to  the  spherical  surface  Sj,  of  the  spherical 
composite  element.  Because  of  the  absence  of  radial  symmetry  the  boundary 
value  problem  to  be  solved  is  much  more  difficult  than  that  of  the  preceding 
section.  Its  formulation  is  as  follows: 
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Solve  the  two  systems  of  field  equations  of  the  theory  of  elasticity: 

(4.1) 


(Xp  +  Gp)e|P^  +Gp  0 


(p)_  „(P) 


6  =  U 


(V  +  V  =  0 


(m)  (m) 


(4.2) 


=  '^i.i 


} 


where:  i  -  1,  2,  3,  a  comma  denotes  partial  differentiation,  a  subscript  or 
superscript  p  denotes  the  inclusion,  a  subscript  or  superscript  m  denotes 


the  shell,  ^ 


The  boundary  conditions  are 


denotes  Laplace' s  operator. 


u<P)  .  0 

1 


,{P)  (m) 

Ui  =  u^ 


Ti 


(P) 


,(m) 


(m) 


=  rn- 


(r  =  0) 


(r  =  a) 


(4.3) 

(4.4) 

(4.5) 


=  rn. 


=  0 


(r  =  b) 


(4.6) 


Equations (4. 6)  express  the  condition  that  pure  shear,  is  applied 

to  the  boundary  in  the  stress  approach.  The  direction  cosines  n^  are 
given,  for  a  spherical  surface,  ly 


“i 


=  5L 


(4.7) 
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In  the  displacement  approach  a  homogeneous  shear  displacement 
is  applied  to  the  boundary.  Then  (4.6)  is  replaced  by: 


(m) 

(m) 

u. 


(r  =  b) 


(4.8) 


An  exact  closed  solution  in  terms  of  solid  spherical  harmonics  is 
given  in  the  appendix  to  this  paper.  Proceeding  in  a  way  analogous  to  that 
of  section  3,  the  foUcwing  expressions  for  the  approximate  bounds  for  the 
shear  modulus  are  found  (see  appendix  A) 


* 

G| 


'm 


1+  (1  -7)yf^c 


G*  =  G^[l  +  (-?  *1)71'^  c] 


(4.9) 

(4.10) 


in  which 


7  =  — 

Gm 


and  y|°^  and  y^^^  are  defined  b7  equations  (A-28)  and  (A-31)  in  appendix  A. 
So  according  to  (2.44) 


of  ^  G  * 


As  for  the  bulk  modulus  the  results  can  be  easily  generalized  to  the 
case  of  different  kinds  of  particles.  Using  the  following  notation: 
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Gp 

'7i 


fVk 

shear  modulus  of  i  kind  of  particles 

.(i) 


'm 


^3 

)  =  c 


i  =  k 

-  Z  e. 


i  =  l 


where  k  as  before,  is  the  number  of  different  kinds  of  particles,  equations 
(4.9)  and  (4.10)  are  now  replaced 


g'  >  Gm 


^2  =  Gm 


1  =  K 

z 


W 


i=i  i  +  - ’^Oyid)  *^1 


r  i=k 

•5". 


1 


(e) 

1)  71  (i)  Ci 


(4.11) 


(4.12) 


The  and  are  determined  from  (A-28)  and  (A-31)  for  the  values 


.(i) 


Gp^  and  V  .  The  value  of  /  is  not  affected  and  is  as  given  above. 

The  results  for  the  shear  modulus  are  much  more  complicated  than 
the  simple  expression  for  the  bulk  modulus  obtained  in  section  3.  Whenever 
the  bounds  for  the  shear  modulus  are  close  together  a  formula  which  gives 
values  which  lie  between  them  can  be  regarded  as  a  good  approximation.  A 
simple  formula  of  this  kind  is  derived  in  appendix  B.  The  result  is: 

^  15(1- Vni)(~  -l)c 

p-*  V  Gjh  / 

Gm 


7-5Vp,  +  2(4-5Vj^) 


92,  . 


*m 


*m 


1  c 


(4.13) 
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When  there  are  k  different  kinds  of  particles 


(4.14) 

It  is  shown  in  appendix  B  that 


cf  6  (4.15) 

In  view  of  the  numerical  calculations  and  comparison  with  experiments 
which  will  be  given  in  the  next  section  it  is  necessary  to  determine  what  boimds 
on  Young' s  modulus  are  given  by  the  bounds  on  the  shear  modulus.  It  is  easily 
proved  that  an  upper  or  lower  botmd  on  the  shear  modulus  gives  the  corresponding 
bound  for  Young  ‘  s  modulus.  Assuming  that  Gj  and  Gg  are  lower  and  upper 
bounds,  i.e. 

G2  >  G,  (4.16) 
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E, 

■*.>u 

Et 


_ 


1  + 


3K 

G, 


>  1 


from  inequaliiy  (4.16).  So 


E,  >  E, 


(4.17) 


It  is  easily  proved,  in  a  similar  way,  from  the  expression  of  Poisson' s 
ratio  the  bulk  and  shear  moduli,  that  the  upper  boimd  to  the  shear  modulus 
corrftspofids  to  a  lower  bound  on  Poisson' s  ratio  and  vice  versa. 


5.  jsxpressions  for  the  moduli  for  very  small  and  very  large 

CONCENTRATIONS 

A  small  concentration  formula  for  a  modulus  of  a  heterogeneous 
material  is  one  that  is  valid  when  squares  and  higher  powers  of  the  volume 
concentration  of  the  pai’ticles  can  be  neglected.  Thus  for  very  small  c 


•■m 


=  1  +  ac 


=  1  +  ^c 


(5.1) 


(5.2) 


The  numbers  a  and  ^  may  be  interpreted  as  the  slopes  of  the  curves 
"modulus-ccncentration"  at  c  =  0 . 

Similarly  large  concentration  formulae  can  be  defined  by 

S-  5=  1  +  a'  c'  (5.3) 

Kp 

1+^'c' 

Gp 


(5.4) 
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in  which 


c'  =  1  -  c 


(5.5) 


and  c  is  very  close  to  unity. 

The  value  of  a  is  easily  foxmd  from  (3.20)  when  the  term  containing 
c,  in  the  denominator,  is  neglected.  The  following  expression  is  obtained. 


For  the  shear  modulus  equations  (A-2G)  and  (A-31)  have  to  be  solved 
for  the  case  where  f^=  c  is  very  small.  The  upper  and  lower  bound 
expressions  for  small  concentration  are  then  found  from  (A-35)  and  (A-36). 
It  is  found  that  both  expressions  are  the  same  and  reduce  to: 

1 

5 


,  15  (1  -Hn)  ( 

a_  =  1  - - 

Gm 


7-5V^  +  2(4- 


Gm  ^ 


(5.7) 


m' 


*m 


Thus  the  upper  and  lower  boimd  have  the  same  slope  at  c  =  0 .  Expressions 
(5.S)  and  (5.7)  have  already  been  given,  in  the  same  form  by  Hashin  (1958)  and 
in  another  form  by  Eshelby  (1957). 

The  same  procedure  can  be  used  for  very  large  concentrations.  The 
results  are: 


=  1  - 


Km'i 

Kp  ) 


2(l-2Un)  +  (l+N^„)-£ 


‘■m 


3  (1  -  ^m) 


(5.8) 
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7-5Vm  +  2(4-5V^) 


15  (1  -  ^4i) 


(5.9) 


Again  the  expresdons  for  the  shear  modulus  are  the  same;  so  the  upper  and 
lower  bound  have  the  same  slope  also  at  c  =  1 . 

It  is  worthwhile  to  note  that  there  is  a  simple  relation  between  a  and 
a'  defined  in  (5.1),  (5.3)  and  ^  and  defined  in  (5.2)  and  (5.4).  From  (5.6) 
and  (5.8) 


From  (5.7)  and  (5.9) 


(5.11) 

4 


Formulae  for  Young' s  modulus  for  small  and  large  concentration  can 
now  be  easily  derived.  The  relation  between  Young' s  modulus  E  and  the  bulk 
and  shear  modulus  is  given  by: 


E  = 

3  K  +  G 


(5.12) 


Introducing  (5.1)  and  (5.2)  into  (5.12)  and  linearizing  with  respect  to  c  gives 

E^.  =  1+  3%^+Gjaa  ^ 

Em  3  K|jj  +  Gjjj 

Similai  iy  for  large  concentration,  using  (5.3)  and  (5.4), 

si  =  1  +  ^  Kp  ■Hip  {5.14) 

Ep  SKp.Cp 

The  quantities  a,  a'  and  /3'  are  given  in  (5.6)  *  (5.9). 


1 


-2'J- 

6.  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS 


The  theoretical  results  obtained  above  will  now  be  used  for  a 
numerical  comparison  with  experimental  measurements. 

A  suitable  example  is  given  by  Tungsten  Carbide -Cobalt  alloys  (WC-Co) 
the  micro-structure  and  mechanical  behavior  of  which  have  been  extensively 
studied.  The  most  recent  experimental  work  on  this  allc^  is  given  in  papers 
by  Gurland  (1959)  and  Nishimatsu  and  Gurlaad  (1959). 

The  alloy  prepared  consisted  of  WC  particles  imbedded  in  a  matrix  of 

V 

Co.  The  experimental  results  of  chief  interest  for  the  present  work  are 
measurements  of  Young' s  modulus  of  the  alloy  for  varying  volume  concentrations 
of  one  phase  relative  to  the  other.  These  results  are  given  in  Table  No.  1. 
Another  quantity  of  interest  is  the  contiguity  which  is  defined  as  the  average 
fraction  of  surface  area  shared  by  a  grain  of  WC  with  all  neighboring  grains 
of  the  same  phase.  It  has  been  assumed  in  the  theoretical  analysis  that  the 
contiguify  is  zero  for  all  concentrations.  This  can  clearly  not  be  expected  to 
hold  for  an  actual  alloy.  Table  No.  1  contains  the  measured  contiguities.  The 
volume  concentration  c  refers  to  the  WC  phase. 

The  values  of  the  contiguities  show  that  there  is  a  definite  preference 
for  the  Co  phase  to  be  classified  as  the  matrix.  This  is  also  shown  by 
photographs  of  the  microstructure  of  the  alloy,  at  different  concentrations, 
which  are  given  in  the  above  cited  papers. 

In  accordance  with  the  preceding  analysis  the  moduli  of  the  WC  will 
be  given  the  sqbscript  p  and  those  of  the  Co  the  subscript  m  .  Then  in  the 
present  case 


Ejjj  =  30  X  10®  psi 
Vqj  =  0.30 


Ed  =  102x10®  psi 


'P 

Vp  = 


i! 


i 


?! 


0.22 


Using  the  formulae 
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K 


£< 

3(1  -2V) 


(6.1) 


G 


E 

2(14  V) 


(6.2) 


the  values  of  the  bulk  and  shear  moduli  for  the  two  phases  are: 


[m  =  25.0  X  10®  psi 

Kp  =  60.7x10®  psi 

rjjj  =  11.5  X  10®  psi 

Gp  =  41.8x10®  psi 

The  variation  with  volume  concentration  of  the  bulk  modulus  K  ,  of 
the  bounds  for  the  shear  modulus  Gj  and  and  of  have  been 
determined  the  method  given  in  the  theoretical  part  below.  Young's 
modulus  and  Poisson' s  ratio  have  then  been  calculated  from  the  expresstons: 


3K4G 


(6.3) 


3K  -  2G 
2  (3K4G) 


(6.4) 


In  (6.S)  and  (6.4),  K*  given  by  (3.20)  is  used  for  K.  The  bounds  on  E* 
and  V  are  obtained  when  introducing  the  bounds  and  G^  for  G . 
The  quantities  E  and  V^,  which  are  obtained  when  using  G  (given  1^ 
(4.13) )  for  G ,  may  be  regarded  as  approximate  values  for  the  Young' s 
modulus  and  the  Poisson' s  ratio  of  the  heterogeneous  material,  whenever 
the  bounds  are  close.  The  results  are  given  in  Table  No.  2  below  in  form 
of  ratios  to  the  matrix  moduli. 
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Figure  4  shows  the  calculated  variation  of  the  bul':  modulus  of  -die 
alloy.  Figure  5  gives  the  boimds  for  the  shear  modulus  and  the  variation 
of  G  ,  Fig.  6  -  the  bounds  for  Poisson' s  ratio  and  V  ,  Fig.  7  -  the 
bounds  for  Young' s  modulus  and  E*'^.  In  the  variations  of  the  moduli 
results  of  formulae  for  small  and  large  concentration,  given  in  Section  5 
are  shown  as  straight  lines  tangent  to  the  curves  at  concentrations  0  and 
1.0. 

It  appears  from  the  figures  that  the  bounds  are  quite  close  together 
and  are  themselves  good  approximations  to  the  values  of  the  moduli.  The 
experimental  results  given  in  Table  1  are  plotted  in  Fig.  7.  This  shows 
that  the  experimental  points  closely  follow  the  curves  predicted  by  theory. 

7.  DISCUSSION  AND  CONCLUSIONS 

It  has  been  shown  that  approximations  to  bounds  for  the  elastic 
moduli  of  composite  materials  can  be  determined  by  use  of  the  variational 
theorems  of  the  theory  of  elasticity.  For  a  specific  example  of  an  alloy  for 
which  experimental  data  are  available  it  was  found  that  the  theoretical 
results  are  close  to  those  found  by  experiment. 

The  mathematical  expressions  for  the  bounds  for  the  shear  modulus 
depend  in  a  rather  complex  maiuier  on  the  ratio  between  the  shear  moduli 
of  the  two  materials  and  their  Poisson' s  ratios.  It  is  therefore  difficult  to 
give  a  general  criterion  for  the  closeness  of  the  bounds. 

A  possible  procedure,  which  is  applicable  to  specific  cases,  is  the 
A  given  by, 
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following:  The  quantity 


=  25.0  X  10  psL 


SMALL  CONCENTRATION  FORMULA 
LARGE  ■■  11 
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VARIATION  OF  BULK  MODULUS  WITH 
CONCENTRATION 


FIG.  5  APPROXIMATE  BOUNDS  FOR  SHEAR  MODULUS 
AND  G* 


APPROXIMATE  BOUNDS  FOR  POISSON'S 
RATIO  AND  V* 


APPROXIMATE  BOUNDS  FOR  YOUNG'S  MODuLU 
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is  evaluated  numerically  at  a  concentration  of  50  per  cent.  This  will  give 
a  good  estimate  of  the  maximum  difference  between  the  bounds.  Whenever 
A  is  small  enough,  G  given  by  (4.13)  can  be  used  as  the  shear  modulus 
of  the  heterogerieous  material. 

It  is  to  be  expected  that  whenever  the  difference  in  moduli  between 
particle  and  matrix  material  is  not  too  large,  the  bounds  will  be  close 
together.  It  shoidd  be  noted  that  for  the  WC-Co  alloy  the  ratio  ^  between 
the  shear  moduli  was  3.62  which  is  certainly  not  small.  The  bounds  were 
nevertheless  close  together.  It  is  consequently  believed  that  the  simple 
expression  G  can  be  used  with  good  accuracy  for  many  practical  cases. 

A  different  situation  may  arise  for  such  extreme  cases  as  rigid  particles 
or  empty  cavities.  This  is  further  discussed  below. 

An  important  property  of  the  approximate  boimds  should  be  emphasized. 
At  both  extremities  of  the  concentration  range,  c  =  0  and  c  =  1  ,  the  slope  is 
the  same  for  both  the  upper  and  lower  bounds  for  any  elastic  modulus  of  the 
heterogeneous  material.  This  provides  some  additional  foundation  to  small 
concentration  theory  and  introduces  large  concentration  formulae  which,  it  is 
believed,  are  of  the  same  order  of  accuracy  as  those  for  small  concentration. 

Some  consequences  of  a  basic  assumption  made  in  this  work  should  be 
pointed  out.  It  has  been  assumed,  that  at  any  concentration,  a  particle  can 
always  be  surrounded  a  surface  which  lies  wholly  in  the  matrix.  In  other 
words,  the  matrix  always  remains  connected.  A  possible  geometrical 
configuration  of  such  a  material  at  high  concentration  is  shown  in  Fig.  8a. 

Another  possibility  is  a  heterogeneous  material  which  is  an 
agglomeration  of  grains  of  two  or  more  different  materials.  A  concentration 
of  100  per  cent  is  reached  when  all  grains  are  replaced  1^  grains  of  one 
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material.  The  geometrical  configuration  of  such  a  material  at  high 
concentration  is  shown  in  Fig.  8b.  In  this  second  case  there  is  no 
preference  of  one  material  over  the  other.  The  expressions  for  the  moduli 
have  to  be  invariant  for  replacement  of  one  material  by  another  and  of  c 
by  1-c. 

The  analysis  given  in  this  paper  applies  only  to  the  first  case  and 
it  should  be  expected  that  for  two  such  materials,  composed  of  the  same 
constituents,  the  difference  in  moduli  would  increase  with  increasing  volume 
concentration.  This  is  illuminated  by  the  fact  that  while  small  concentration 
formulae  for  these  two  materials  are  the  same,  the  large  concentration 
formulae  are  not. 

The  two  materials  described  should  be  regarded  as  extreme 
theoretical  cases.  The  structure  of  an  actual  heterogeneous  material  will 
in  all  probability  be  somewhere  between  the  two. 

The  theory  developed  in  this  work  may  be  easily  applied  to  the  special 
cases  of  a  body  containing  empty  cavities  or  a  body  containing  rigid  particles. 
In  the  first  case 

Kp  =  0  ;  Gp  =  0 

and  in  the  second 


Equations  (A-28)  and  (A-31)  are  greatly  simplified  in  these  two  cases  and  the 
boimds  can  be  easily  expressed  in  closed  form.  Numerical  calculations  show 
that  the  bounds  so  obtained  are  further  apart  then  for  the  alloy  treated  in 
this  work. 
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The  case  of  cavities  is  important  for  the  study  of  porous  materials. 
The  solution  for  rigid  particles  can  Iqr  mathematical  analogy  be  applied  to 
the  problem  of  the  viscosity  of  a  suspension.  As  there  is  complete 
ma^Jiematical  analogy  between  the  theory  of  incompressible  elastic  media 
ar.d  the  theory  of  Stokes  flow  (viscous  flow  when  neglecting  inertia  terms  in 
the  Na-vier -Stokes  equations),  the  results  for  the  shear  modulus  of  an 
incompressible  elastic  body  containing  rigid  particles  will  hold  for  the 
coefficient  of  visco'sity  of  a  viscous  fluid  in  slow  motion,  containing  rigid 
particles.  However  for  both  cavities  and  rigid  particles  a  special  difficulty 
appears  -  the  range  of  volume  concentration  will  not  extend  to  100  per  cent. 
When  the  skeleton  of  a  porous  material  breaks  down  the  volume  of  solid  oi 
which  this  skeleton  is  composed  is  not  negligible  in  comparison  to  the  voliune 
of  the  voids.  Also  the  flow  of  a  suspension  will  cease,  i.e.  -  the  coefficient 
of  viscosity  becomes  mfinite,  at  a  concentration  of  50  per  cent  -60  per  cent 
of  partick's. 

It  seems  that  the  method  given  in  this  work  cannot  be  applied  without 
modification  to  these  extreme  cases. 
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APPENDIX  A 

PROBLEM  OF  SHEARED  COMPOSITE  SPHERE 
AND  BOUNDS  FOR  SHEAR  MODULUS 

In  the  following  a  solution  to  the  boundary  value  problems,  set  up  in 
section  4,  will  be  constructed. 

General  solutions  for  the  elastic  sphere  and  the  elastic  spherical 
shell  have  been  given  by  Lord  Kelvin  (Thomson)  and  Tait  (1879),  and  Lure 
(1953)  for  the  general  case  and  by  Sternberg,  Eubanks  and  Sadowsky  (1957)  for 
the  axisymmetric  case .  It  is  in  principle  possible  to  construct  solutions  for  the 
spherical  inclusion  and  for  the  spherical  shell,  leaving  enough  arbitrary  constants 
to  fit  the  boundary  conditions. 

The  first  two  general  solutions  are  very  difficult  to  use  owing  to  uie 
complexity  of  the  general  solution  for  the  spherical  shell.  The  axisymmetric 
solution  is  more  explicit  and  it  can  be  applied  to  the  present  problem  by 
superposition  of  axisymmetric  tension  and  compression  at  right  angles.  A 
simpler  method  of  solution  will  be  found  by  proceeding  in  a  less  straightforward 
way. 

The  problem  of  a  spherical  inclusion  in  an  infinite  elastic  mediiim,  the 
moduli  of  which  are  different  from  those  of  the  inclusion,  when  the  stress  at 
infii\ity  is  imiform  (or  the  displacements  are  linear)  may  be  solved  in  a 
convenient  way  by  application  of  Kelvin' s  general  solution  of  the  problem  of 
an  elastic  sphere  [  (Love)  pp.  265-270] . 

In  the  case  of  pure  shear  in  the  1,  2  directions  it  is  found  for  the  region 
exterior  to  the  inclusion  that  displacements  and  accordingly  stresses,  may  be 
expressed  in  terms  of  a  solid  spherical  harmonic  of  negative  integral  degree 
*3, 
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n-3)  -  ^ 

and  its  derivatives. 

Love  (pp.  249-251)  has  given  three  special  solutions  to  the  field 
equations  of  elasticity  and  it  can  be  shown  that  the  solution  of  the  problem  just 
mentioned  may  also  be  found  superposing  Love's  w  and  ^  solutions. 
Ttiese  are  in  vector  form: 

Type  tt) 


tT  =  r^  grad  +  %  “n  ^ 

an  =  -  2  nX-H(3n-t-l)G 
(n  +  3)  X  +  (n  +  5)  G 

Type 


(A-2) 


vT  =  grad 

For  the  infinite  medium  in  pure  shear  the  solution  can  be  constructed 
by  putting 


“n  =  '#'n  -  -3) 

In  the  problem  given  by  equations  (4.1)  -(4.8)  the  exterior  region  is  a  finite 
spherical  shell.  It  is  reasonable  to  expect  that  the  solution  can  be  constructed 
by  using  also  the  solid  harmonic  of  positive  integral  degree  associated  with 

^(-3)' 

If  is  a  solid  harmonic  of  positive  integral  degree  n  then: 


j.2n+l 


'^^-n-1  = 
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is  a  solid  harmonic  of  integral  degree  -n-1.  [Compare  e.g. 

MacRobert 

(1948)  p. 

74.]  So  in  this  case 

^3  =  X2 

(A-4) 

The  elastic  displacement  in  the  spherical  shell  can  now  be  written  as 

follows, 

=  A^TT^^^  +  A2 

(A-5) 

where, 

=  grad  (Xj  X2) 

(A-6) 

T?<2)  =  grad("^"^) 

(A-'O 

grad  (Xj  X2)  +  <*1™^  Xj  X2 

(A-8) 

-»(4)  2  j  Z'  *2  ^  *2  ■z^ 

^  ♦  “(-3)  ^ 

(A-9) 

According  to  (A-2) 

(m)  2Xni  +  7  Gm 

a~  -  -  2 

5X^+  7G^ 

(A-10) 

(m)  3  Xjji  +  8  Gjjj 

®/  .X  - 

(A-11) 

Simiiarly  for  the  interior  of  the  inclusion 

“u  =  B,  +  B3  (A-12) 

Bj  =  Bif  =  0 


(A-13) 
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The  last  equation  is  derived  from  the  fact  that  u' 
For  the  interior  region 

^(3)  1.2  ^^(p)  ^ 


(P) 


=-2^ 


2  Xn  +  7  Gr 


5  Xp  +  7  Gp 


vanishes  for  r  =  0 . 


(A-14) 

(A-15) 


In  the  fcUov^ing  expressions  for  the  stress  vectors  on  a  spherical 
surface  will  be  needed.  These  may  be  determined  by  differentiation  of  the 
displacements  and  are  given  for  the  special  type  solutions  in  Love 
(pp-  250-251).  Inti'oducing  the  nondimensional  constants 


Ax 


A^ 


A3  = 


A 


f 


(A-16) 


Bx  =  B, 

Bj  =  a  Bj 


and  using  the  formula 
A  =  2^ 

G  1  -2V 

where  V  is  Poisson's  ratio  -  the  stress  vectors  on  a  sphere  of  radius  r 
where  a  <:  r  5;  b  are  given  by  the  following  expressions: 
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Ml  1-2V^ 


7+5V 

+  20  A,  -  2  - P- 

'  7-4V^ 


m 


(A-17) 


Tf>=  IPm  A  .4^  !l22k  Aj+  A, I  X, 

^  L  1.2V^ 


’7+5Vm  12  "  X,  xi 


(A-18) 


(A-19) 


The  remaining  components  are  easily  written  down  by  analogy  with  (A-18)  - 
(A-19). 

The  expressions  for  the  displacements  are;  For  the  shell 
a  ^  r  ^  b 


-38 


=  (Aj  +  Ag  +  Aj  +  A^)  X2  + 

./  7-10V„,  .  3 


;  I  -AU  y'rn  *<  \  X* 


2 

Xf  X2 


(A-21) 


(Aj  +  A2  +  As  +  Ai^)  Xj  + 


(- 


+  I  -  5  A2  -  2  - 


7  “  10  V, 

-4Vjjj  1  -2V~^V 


Iffi  A  .  _ 3  _  ^  ^  Xa 


(A-22) 


u<f>  =  (-5A,.2P^A..^ 


V  ‘*"2  ^  7  -  4  Vjjj  1  -  2  Vjjj  "“t 


X2  X3 


(A-23) 


For  the  interior  of  the  inclusion,  0  &  r  i  a  : 

,(P).,„  B, 

T*  4V  r 


u\  =  (B1+B3)  Xg  -  2 


(A-24) 


and  the  other  components  follow  by  analogy. 

The  two  boundary  value  problems  for  the  elastic  composite  sphere 
may  now  be  easily  solved.  The  expressions  given  for  the  displacements 
already  satisfy  equations  {4.1)j  (4.2).  Equation  (4.3)  is  satisfied  by  the 
components  of  u^^  given  by  (A-24).  Setting  {A-21)  equal  to  (A-24)  in 
accordance  with  (4.4),  (A -17)  equal  to  (A -20)  in  accordance  with  (4.5)  and 
equating  to  zero  coefficients  of  Xg  and  x/  Xg  -  4  linear  equations  for  the 
six  unknown  coefficients  A,,  Ag,  Aj,  A^,  B^  ,  B3  are  obtained. 

It  should  be  noted  that  eoiiating  of  other  components  of  displacement 
and  stress  vectors  does  not  give  new  equations.  This  follows  from  the  form 
of  the  expressions  and  may  easily  be  understood  by  the  symmetry  of  the 
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problem  with  regard  to  the  bisector  plane  drawn  between  the  Xj  Xj  and 
Xg  Xj  planes. 

Two  additional  equations  are  obtained  for  each  of  the  problems,  either 
from  (4.6)  and  (A-17)  or  from  (4.8)  and  (A-21).  The  problem  is  thus  reduced 
to  the  solution  of  six  linear  equations  with  six  unknowns. 


For  the  stress  problem  the  matrix  of  coefficients  is  as  follows: 


and  f  has  been  defined  by  (3.1)  and  (3,2). 
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For  the  displacement  problem  only  rows  (5)  and  (6)  in  (A-25)  undergo 
a  change.  These  are  now: 


The  system  of  equations  {A-25)  may  be  reduced  after  various 
rearrangements  and  application  of  rules  for  evaluation  of  determinants  to 
a  system  of  two  equations  in  two  unknowns: 

yP  1 1  if'  +y^?|[(7-10Vp)-(7-10V^)a^4/-  (7+5x41)^^!  =  0 


yx  n 


15  <1  -Vj„) 


(1-7)  (1-7) 


+  72  ^  i  j 


(A-28) 


in  which. 


>)  =  ISSL 

71  ^ 


21 


5(7- 4  Vp) 


(0)  (o) 

yy  =  sr 


(A-29) 


_  ^(7-  10^p)+7  (7-»-5Vp) 
“  35  (1  -  VJ 


(A-30) 


V  .  S’ 


\ 


The  system  of  equations  associated  with  (A*27)  may  be  reduced  in  a  similar 
way  to: 

y?  •  I  (/-/)}  +y?  [(7*10Vp)-(7-10Vjjj)-l^4/+4{7-10^^jn)^=  0 

(A-31) 


in  which 


I  \^K 


5(7-4^Jp) 


(A-32) 


and  is  given  by  (A -30). 

The  quantities  i  S  are  best  determined  in  this  case  from 
equations  (2.14)  (2.13)  using  (A-2I)  and  Us  analogues,  and  (A-24)  and  its 
analogues.  The  unknown  constants  B,  ,  Bj  are  determined  from  (A-28) 
and  (A -29)  for  the  stress  approach  and  from  (A-31)  and  (A-32)  for  the 
displacement  approach.  The  ^o) 

and  uj°^  appearing  in  (2.14)  are  given 
by  (4.6)  and  (4.8)  respectively.  (Because  of  the  theorems  that  when  an 
elastic  body  is  loaded  by  constant  stress  this  stress  system  is  found  at 
every  point  in  it;  when  deformed  by  linear  surface  displacements  this 
displacement  field  is  found  throu^out  it.)  Carrying  out  the  integration  the 
following  expressions  are  foimd: 


7 


\ 


\ 


W'*  =^(i->!)yf  V„ 


(A-33) 


..(e)  ^ 


^  - 1)  y\'  Vn 


(A-34) 


Then  from  (2.45)  and  (2.46) 


G,  = 


1  +  (1  -  -7 ) 


(A-35) 


=  Gm[l  +  (^  -l)y?  c] 


(A-36) 


where  (2.47)  and  (2.22)  have  been  used. 


The  actual  expressions  for  and  are  clum^  and  it  seems 


best  to  determine  them  numerically  when  numerical  values  of  Kjj^,  Gm 

.  ...  -  -  _ . _  ...  _ U  1 _ _ J  «^^)  n(  )  n  wn  fitvtni 


and  Kp,  Gp  are  given.  It  should  be  noted  that  y\  and  y;'  are  functions 
of  f  i.e.  of  the  voliime  concentration  c. 


V..'  — 'r'x 


A 


*7 


.APPENDIX  B 

APPROXIMATE  EXPRESSION  FOR  SHEAR  MODULUS 


The  expression  (4,13)  for  G  can  be  derived  in  the  following  way. 
The  system  of  equations  (A-28)  is  written  in  the  form 


(a)  (j)  (0) 

ail  yi  +  ai2  ya  =  0 


^  e  IP  -  1 


(B-1) 


And  analogously  for  system  (A -37), 

(c)  (e)  (6)  (€)  . 

^11  yi  +  ai2  ya  =0 

(e)  (e)  (e)  (e)  , 

^21  Yi  +  y^  =  1 


(B-2) 


Then, 


^11  ^22 


V/ 


(«) 

Yi  = 


(^)  ajj  ags 


CB-4) 


The  quantities  and  y|^^  are  defined  ly, 

-(c^)  _L_ 

=  >) 


(B-5) 


03-6) 
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If  and  yj'  in  equations  (4.9)  and  (4.10)  [  (A-35)  and  (A-36) 
in  Appendix  AJ  are  replaced  by  and  yf^^ ,  the  resulting  expressions 
are; 


gT  ,  1 

1  +  (1 -'’^)  y^'^^  c 


(B-7) 


£2_  =  l  +  ('>^  -1)  yj^^  c  (B-8) 

Gm 

Introducing  the  coefficients  and  from  (A-28)  and  (A-31)  it 
turns  out  that, 

Gr=G3*  =  G'"  (B-9) 

where  G*  is  given  by  equation  (4.13), 

It  can  be  proved  that  for  all  values  of  c, 

Gj  6  G’*'  ^  Gg  (3*10) 


The  proof  is  somewhat  tedious  and  only  its  outline  will  be  given  here. 
This  is  as  follows: 

(a)  The  quantities  y^j  ^  and  y^f'  are  always  positive. 

From  the  definition  of  6  and  S  it  follows  that  when  the 
moduli  of  an  inclusion  are  larger  than  those  of  the  mediiun  <5"  is  negative 
and  <5  is  positive.  When  the  moduli  of  an  inclusion  are  smaller  than  those 
of  the  medium  the  converse-  is  true.  Applying  this  to  equations  (A-33)  and 
(A“34),  statement  (a)  follows. 
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(b)  The  following  relations  hold  for  the  coefficients  and  . 


H  (€) 

^11  "  ^11  "  ^11 

(<r)  (e) 


^22  ^22  "  ^22 


^21  "*  ^21  ~  (1  *7)  C 


(e)  ('3’) 

a/a  -  a,2  =  35(1  -V^)T>' 


(B-11) 


aii  >  0 
a„  <  0 
i<r) 

ax2  ^  0 
(<3) 

aJ  >0 


as  2  ^  ® 


ai2  >  0 

(€) 

^21 


for  y  >  1 
for  ->7  <  1 


V 


>  tB-12) 


(c)  The  quantities  and  satisfy  the  following  inequalities; 

Tf  >  yf 


for 

1 

4  M\ 

for 

-q  <  1 

(B-14) 

This  follows  from  (B*3)  and  ^*4)  and  the  inequalities  given  in  (b). 


(d)  The  approximate  bounds  satisfy  the  necessary  condition: 

of  ^  G3*  (B-15) 

This  follows  from  (A -35),  (A-36),  (B-3),  (B-4),  (B-12),  (B-13)  and 

{B-14). 

(e)  The  quantities  and  G*  are  related  to  and  Gg  by  the  following 
inequalities: 

G*  £  G,^  (B-16) 

G*  6  G*  (B-17) 


This  follows  by  introducing  (B-5)  and  (B-6)  into  (B-7)  and  (B-8), 
(B*3)  and  (B-4)  into  (A-35)  and  (A-36)  and  then  comparing  g/*  with  Gj* 
and  Gg  with  G^*,  using  (B-12). 

So  according  to  (B-16),  (B-17)  and  (B-9)-  (B-10)  is  true. 

The  proof  remains  essentially  unchanged  for  the  case  when  the 


particles  are  of  different  kinds. 

Equation  (413)  may  be  rearranged  into  a  symmetric  form  analogous 
to  (3.21), 

K'-tl 


/  r*  \  ^ 


C  = 


lOL 

N 


(B-18) 


/  G"\/  Gjj  ) 

1  g  +  pr-  ! !  1  - 
V.  ^m/  \  ''m 


In  which 


g  = 


7  -5'' 


m 


2(4 -5Vj^) 

Note  also  that  (4.13)  satisfies  the  slope  conditions  (5.7)  and  (5.9), 


as  it  should. 
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